Two-photon quantum walks in an elliptical direct-write waveguide array 
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Integrated optics provides an ideal test bed for the emulation of quantum systems via continuous- 
time quantum walks. Here we study the evolution of two-photon states in an elliptic array of 
waveguides. We characterise the photonic chip via coherent-light tomography and use the results 
to predict distinct differences between temporally indistinguishable and distinguishable two-photon 
inputs which we then compare with experimental observations. Our work highlights the feasibility 
for emulation of coherent quantum phenomena in three-dimensional waveguide structures. 



I. INTRODUCTION 

Computer simulations of complex systems have had 
spectacular success in modern science due to sophis- 
ticated classical approximation methods and steadily 
increasing computational power. However, it is ulti- 
mately impractical to accurately simulate the dynam- 
ics of even moderately-sized quantum systems due to 
their exponentially- increasing parameter space. Quan- 
tum simulation, as first proposed by Feynman [1 , offers 
a path forward whereby the simulator itself operates via 
quantum instead of classical dynamics. 

There are two possible approaches to modeling a quan- 
tum system. The first is simulation^ in which a digi- 
tal model of a system's Hamiltonian yields a physical 
quantity. An example is provided by recent quantum 
chemistry simulations of a Hydrogen molecule [2 [3] on 
a small-scale quantum computer. The second is emula- 
tion^ which involves the modeling of a quantum system 
with a better-controllable system which is governed by a 
sufficiently similar Hamiltonian [4 . 

Quantum walks — natural extensions of the ubiquitous 
classical random walk into the quantum world O |6] — 
provide an ideal framework for emulation due to their rich 
dynamics There are two limiting cases, discrete 

and continuous. In discrete-time quantum walks, one 
or more (interacting) quantum particles (the so-called 
quantum coin) evolve on a graph, with their evolution 
governed by their internal quantum (coin) states. The 
discrete-time quantum walk on a line is the best studied 
example of such a walk and it has been demonstrated in 
a number of physical systems p!QHT4] . 

In continuous-time quantum walks, in contrast, there 
are no coin operations, the evolution is defined entirely in 
position space [15 . These walks require a well-controlled, 
continuous coupling between vertices, or lattice sites. 
Ideally, this coupling should be stable and robust against 
decoherence. Integrated optics is perfectly suited for this 



task. Lithographically written, evanescently coupled sur- 
face waveguides have already been used to study contin- 
uous time quantum walks on a line [16 . A disadvantage 
of surface waveguides is that they can only realise simple, 
one-dimensional graphs with limited interconnectivity. 

Physically more interesting three-dimensional struc- 
tures can be engineered in laser-written optical waveg- 
uide arrays in dielectric materials such as fused silica [17J. 
The two-dimensional graphs that can be realised with 
this technique allow the study of new quantum walk 
phenomena, such as wave communication [18], cooper- 
ative quantum games and the creation of topologi- 
cal phases in two dimensions [20]. Examples of waveg- 
uide array structures which have been manufactured with 
this technique include rings, regular and hexagonal lat- 
tices, X-shapes and triangular shapes [2TT[23 . To date, 
however, these have only been explored with coherent 
light. True quantum effects, and associated genuine 
improvements in efficiency of emulating quantum sys- 
tems will only emerge for non-classical input states [24] , 
such as recently demonstrated experimentally for a one- 
dimensional waveguide lattice in [25 . 

Here, we study mult i- walker continuous- time quantum 
walks in an optical chip containing an elliptical arrange- 
ment of coupled direct-write waveguides. We characterise 
the optical chip via coherent light tomography, effectively 
realising single particle walks and use the results to pre- 
dict correlations in non-classical two-photon walks which 
we compare with experimental observations. Our work 
is an important step towards the emulation of quantum 
systems in three-dimensional integrated photonic archi- 
tectures. 



II. DEVICE DESCRIPTION 

The circuit for the continuous time quantum walk, 
shown in figure [l] consists of six waveguides written into 
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FIG. 1: Schematic of the integrated waveguide circuit and 
associated output pattern. The fan-in region here is scaled 
by a factor of ^ 54 relative to the interaction region along 
the beam direction. At the input the waveguides are equally 
spaced at 127 /im. They then converge via a two-stage fan-in 
to their eventual elliptical configuration, which is shown in the 
CCD camera picture in the inset: the waveguides are arranged 
with equal angular spacing around an ellipse with semi-major 
and semi-minor radii of 10.2 and 7.0 /xm respectively. 



a chip of high-purity fused silica using an ultrafast direct- 
write technique |T7llMll2l]- In this technique, femtosec- 
ond Ti: Sapphire laser pulses tightly focused inside the 
sample yield localized refractive index modifications. By 
translating the sample in three dimensions these index 
modifications can manifest as waveguides tracing three- 
dimensional curves, a process that cannot be replicated 
in conventional lithographic waveguide fabrication tech- 
niques. Our chip was written with a 1 kHz repetition 
rate, 800 nm, 120 fs laser, passed through a 520 /im slit 
and focused with a 40x 0.6 NA microscope objective. 

At the circuit input, the six waveguides are arranged 
along a line with an equal spacing of 127 /im. This allows 
each waveguide to be addressed individually and simul- 
taneously. The waveguides then converge via a two-stage 
fan-in to their final configuration as shown in figure [l] In 
the primary fan-in stage, which occupies the first 8.5 mm 
of the chip, the waveguides follow S-bend curves from 
their linear input arrangement to an elliptical configura- 
tion twice as large in radius as their final configuration. 
In the second fan-in stage, comprising the next 1 mm, fur- 
ther S-bends shrink this array to its final geometry of an 
ellipse with a semi-major axis of 10.2 jam and semi-minor 
axis of 7.0 /im. All the S-bends are of the 'raised-sine' 
form which has been shown to minimise bend loss [28^, 
while the two-stage fan-in configuration was designed to 
minimise coupling (in particular, asymmetrical coupling) 
between non- adjacent waveguides before they reach their 
final configuration. 

The light evolution in such a waveguide array is gov- 
erned by the evanescent inter-waveguide coupling which 
drops off exponentially as a function of the waveguide dis- 
tance. It is theoretically described by a coupled-oscillator 
Hamiltonian, which is described in the Appendix. 



III. OPTICAL CHIP CHARACTERISATION 

The experimental setup is shown in figure |2j Light was 
coupled into the chip via a V-groove array, which houses 
six single-mode input fibres on a line, matching the input 
spacing of the circuit waveguides. We first illuminated in- 
dividual waveguides in the chip using coherent light from 
an 820 nm laser diode, see figure |2|l). The output inten- 
sity profiles were processed in Matlab and converted into 
probability distributions. When compared to the numeri- 
cal simulation obtained with the software suite which was 
used to design the chip, the predicted distribution at the 
circuit output was starkly different to the observed dis- 
tribution as shown in figure [Sj This behaviour prompted 
us to empirically determine the full optical response of 
the circuit using six-port, polarisation- sensitive coherent 
light tomography. 

For this tomography the input polarisation was set 
using a bare reference fibre on top of the chip. Into 
each waveguide we input the following set of polari- 
sation states: {\H) , \V) , \D) , \A) , \L) , \R)}, where \H) 
and \ V) represent horizontal and vertical polarisation and 
\D/A) =m ± \V))/V2, \L/R) ={\H) ± i \V))/^. For 
each output we measured the same six polarisation com- 
ponents, obtaining a total of 216 CCD images. 

Using output intensity distributions taken from these 
images, we subsequently reconstructed an array of 36 
Mueller matrices M [29 . This array completely charac- 
terises the optical response of the circuit, quantifying the 
effects of three distinct processes: notably polarisation- 
dependent inter-waveguide coupling, birefringence, and 
polarisation-dependent loss [29 . 

The results indicate strong birefringence in each of the 
six waveguides. Most notably, when the state \H) is in- 
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FIG. 2: Experimental setup. The chip can be addressed with 
6 individual single mode fibres. The input light polarisation is 
set by a combination of half- and quarter- wave plates (HWP, 
QWP) and polarising beamsplitters (PBS). The chip output 
is magnified and then collimated with two spherical lenses, 
(a) The chip was characterised with an 820 nm laser diode, 
imaged onto a CCD camera via a polarising prism, (b) Quan- 
tum walks were performed with two-photon inputs created via 
spontaneous parametric downconversion (SPDC). The rela- 
tive delay between the two input photons was adjusted 
using a translation stage. We used a 50/50 beamsplitter (BS) 
and two ^ 500 /xm apertures to select a combination of output 
ports and measure the two-photon correlations in coincidence 
using single-photon avalanche diodes. 
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FIG. 3: Comparison of numerical simulations and observed 
probability distributions for the optical chip in figure [l] (a) 
The propagation dynamics predicted using optical waveg- 
uide simulation software (see Appendix) with light input into 
waveguide 1, as a function of z. The simulation includes mod- 
elling of inter-waveguide coupling during the 1 mm second 
fan-in stage at z<0. The curves for waveguides 2 and 4, and 
3 and 5 overlap due to symmetry, (b) The predicted distri- 
bution at the output of the circuit, and the corresponding 
observed probabilities. The distinct differences in the two 
distributions indicate that the optical response of the chip is 
not scalar, instead suggesting significant vectorial behaviour. 
This was confirmed by further tomographic analysis. 



put into waveguide 5, 29% of the overall output state 
across the six channels is rotated to |V). For the input 
channels 1, 2 and 4, the overall polarisation rotation was 
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FIG. 4: Poincare sphere representation of the correspond- 
ing Mueller matrices, Mout,in (see Appendix) describing the 
transformation from input waveguide 6 to outputs 1 — 6. Input 
states shown on the outer Poincare spheres are mapped to dif- 
ferent locations on the blue ellipsoids, simultaneously showing 
polarisation-dependent coupling and birefringent effects in the 
respective channel. The relative orientation of the ellipsoids 
is indicated by the set of orthogonal states \H) (red arrow), 
\D) (green arrow) and \R) (black arrow), and the point of 
contact between sphere and ellipsoid is indicated by the yel- 
low arrow. Note that the arrow lengths are proportional to 
output power, not degree of polarisation. The numbers above 
each sphere give the normalised average power coupled into 
the respective channel. 



small, with more than 91% of \H) being retained in each 
case. Furthermore, there was significant polarisation- 
dependent coupling between the waveguides for all input 
channels. For instance, for input \H) into waveguide 1, 
80% of the total output intensity was observed in channel 
6, however when the input state was |V), only 11% of the 
total output intensity was contained in this channel. An 
exemplary selection of Mueller matrices, illustrated on 
Poincare spheres and quantifying these effects is shown 
in figure [4] The full matrix array A4 can be found in the 
Appendix. 

In addition, the whole chip exhibited significant polar- 
isation dependent loss; integrating over all output chan- 
nels, we observed an excess 38% loss of \V) compared to 
\H) for waveguide 6. This may be due to a combination 
of absorption into the bulk of the circuit, or polarisa- 
tion dependence of the input coupling efficiency from the 
V-groove array to the chip, or both. 



IV. TWO-PHOTON WALKS 

The coherent-light tomography encompasses all possi- 
ble single-photon walks in this chip, since a single photon 
shares the coherence properties of a coherent light beam. 
However, these walks can be efficiently simulated classi- 
cally and it is not until we input multi-photon states that 
we observe quantum effects [24 . 

A schematic of the setup for the two-photon walk ex- 
periment is shown in figure |2]3). Pairs of single photons 
are created via spontaneous parametric down- conversion: 
a mode-locked 76 MHz laser centred at 820 nm is fre- 
quency doubled to 410 nm and focused into a 2 mm long 
/3-barium borate (BBC) crystal, phase-matched for type- 
I downconversion. After passing through interference fil- 
ters at 820 ±2.5 nm, the degenerate output photon pairs 
are coupled into single-mode fibres equipped with man- 
ual polarisation controllers. They are then connected to 
the V-groove fibre array and thence to the chip. 

We created photon pairs at a measured rate of 180 kHz, 
of which we detected a total of 6.5 kHz at the circuit out- 
put. The main loss contributions stemmed from the poor 
input coupling efficiency between the V-groove array and 
the chip (-31%), and the 50/50 beamsplitter (- 50%) at 
the output. The imperfect coupling is mostly due to a 
slight mismatch between the spacing of the fibres in the 
V-groove array and the locations of the input ports of the 
circuit. The observed loss could be significantly reduced 
by using a more sophisticated imaging system. 

We carried out two-photon quantum walks in two sep- 
arate scenarios: with temporally distinguishable and in- 
distinguishable photons. When the photons entering 
the chip are temporally distinguishable, i.e. with a 
time delay larger than their respective coherence times, 
they perform independent quantum walks with local 
evolution. When they enter the chip simultaneously, 
Az = 0, they experience non-classical two-photon inter- 
ference [30 , which can yield quantum dynamics, includ- 
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FIG. 5: Example of observed two-photon interference between 
output waveguides 2 and 6 as a function of relative path dif- 
ference between photon pairs input into neighbouring waveg- 
uides 1 and 2. The visibility of the dip is V2,6 = 38 ± 2%, 
calculated from a Gaussian fit (blue line). 



ing generation of two-photon entanglement. The the- 
oretic qualitative difference between these two cases is 
described in the Appendix. Figure [5] shows an exem- 
plary calibration scan of coincidence counts as a function 
of time delay. The signature of indistinguishable quan- 
tum walkers manifests as a dip in the rate of coincident 
detection events, C, at zero delay, with an interference 

visibility of V = (C^ax - C^in)/C^ax of 38 ± 2%. 

The results for two-photon quantum walks for distin- 
guishable and indistinguishable photons input into the 
neare^^-neighbour channels 1 and 2 are shown in fig- 
ure [6^), as the normalised coincidence probability dis- 
tributions, F^ and F* respectively. Distinct differences 
are observed between the two cases, as suggested by the 
strong two-photon interference signature in figure [5| 

The measured distributions are compared with pre- 
dictions, figure [gJd), which are based on determining 
the components of the waveguide array unitary [/, for 
a particular input polarisation, see Appendix for de- 
tails. The generalised overlap fidelities defined in the 
Appendix, between our measurements and predictions 
are 5*^=93.4 ± 0.3% for the distinguishable walkers, and 
5* =91. 6 ± 0.4% for the indistinguishable walkers. 

Figure [T^) shows correlation matrices for inputs 2 and 
4 as an example of two-photon walks with next-neaiest 
neighbour input ports. We again observe non-classical in- 
terference signatures, with visibilities up to V2,4=28±3%. 
The fidelities between the measured two-photon matri- 
ces and their corresponding predictions, figure [7]3), are 
S'^ = 97.9 ± 1%, and S' = 96.2 ± 0.8%. 



V. DISCUSSION 

In conclusion, we have performed the first quantum 
walk on a general two-dimensional waveguide structure 
with genuine non-classical inputs. This is a significant 
step towards emulating Hamiltonians which can be ap- 
proximated with evanescent ly coupled waveguides. How- 
ever, we have also identified a number of obstacles which 
must be addressed before such an approach can be fully 
realised. 

First, despite the apparent good agreement between 
our predictions and the observed two-photon probability 
distributions, the two-photon visibilities — which quantify 
the measured non-classical effects — matched the predic- 
tions poorly. This is because the Mueller-matrix array 
derived from the coherent beam intensities only yields the 
squared absolute values of the elements of the system's 
unitary U ; it does not determine the (generally complex) 




FIG. 6: Correlation matrices for nearest-neighbour input 
channels 1 and 2. We recorded the photon-coincidence counts 
at each of the 36 pairs of output channels in a 20 second 
time-window, (a) The measured and (b) predicted corre- 
lation matrices for (left) temporally distinguishable photon 
pairs F"^, (center) temporally indistinguishable simultaneous 
walkers F*, and (right) the difference F^ — F\ The coincidence 
probability at the outputs 2 and 6 between the two plots re- 
flects the two-photon interference dip shown in figure [5] The 
measured uncertainties are not plotted, as they are too small 
to be seen on the plots. 
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FIG. 7: Correlation matrices for next-nearest neighbour input 
channels 2 and 4. The plots follow the same convention as 
those in figure [g] for both (a) the measured and (b) predicted 
correlation matrices. Uncertainties are again too small to be 
seen on the graphs. 



phase relations of the waveguide array. These phase fac- 
tors could in principle be obtained by phase-sensitive co- 
herent light tomography, as proposed in [23l [31] . An al- 
ternative technique [32] requires a single A^-photon input 
state (in this case TV = 6) and photon- number resolved 
detection at each output. However, generating large TV 
number states is technically difficult, and both techniques 
suffer from the large number of measurements required 
to fully characterise the six-port system. This is a peren- 
nial problem in quantum science: the exponential power 
granted by multi-photon walks on big lattices makes it 
hard to experimentally characterise system dynamics. A 
potential solution might be to use compressive sensing 
techniques, which have recently been exploited for ex- 
ponentially efficient quantum process tomography [3T. 
Another open question is how to deliberately control the 
phases in the manufacturing step of a waveguide array. 

Second, while we observed significant two-photon in- 
terference visibilities, the resulting probability distribu- 
tions did not exhibit a conclusively quantum signature, 
as quantified by the witness defined in [T6|. This was 
most likely due to the significant polarisation-dependent 
coupling and loss in the circuit, leading to non-unitary 
evolution which failed to preserve the coherence of the 



input quantum state. Once we better understand the 
origins and mechanisms of these effects, we can exploit 
them to engineer devices such as polarisation-dependent 
couplers. It should be noted that these effects are cer- 
tainly not inherent to the direct-write technique, as for 
example demonstrated by [34] . 

Future research should also focus on realising and 
studying effects of decoherence in continuous-time walks. 
Many physical, biological or chemical systems are 
strongly coupled to their environment and decoherence 
has been suggested to be the driving factor behind phe- 
nomena such as environmentally-assisted quantum trans- 
port |9j . While decoherence has been studied in discrete- 
time experiments [14], techniques for controllably intro- 
ducing it to inherently robust waveguide lattices have yet 
to be demonstrated. 
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Appendix 

In the Heisenberg picture, a light field input into a 
waveguide in this chip is subject to the coupled-oscillator 
Hamiltonian 



H 
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where Pi is the propagation constant in waveguide i and 
Cij is the coupling constant between waveguides i and 
j. The system then evolves in time according the unitary 
operator U{t) = e~*^^/^ and the creation operators a] 
are subject to the Heisenberg equation of motion 



dal 
S7 



n da] 



a, 



(2) 



which has the solution 



(e-^K,a](0) = ^f/.,,(z)a](0), (3) 



where C = {Cij} is the 6x6 matrix of coupling constants 
with diagonal entries Ci^i = /3, and z is the propagation 
distance along the waveguide array. Note that this evo- 
lution is equivalent to the continuous-time quantum walk 
formalism [15]. 
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The interaction length z of the waveguides is chosen 
to match the desired run time i for the emulation of the 
Hamiltonian. The overall response of the circuit as a 
12-port beam-splitting device is then contained in the 
unitary matrix U = e*^^, and we can define a set of six 
output creation operators with = aj(z) = ^ UijOj 

3 

for the input operators aj = aj{0). 

The coupling Cij between two waveguides falls off ex- 
ponentially with waveguide separation ^21j, so to de- 
sign an array of the type discussed in this paper, the num- 
ber of waveguides, their geometry, and their separations 
Tij are chosen to refiect the properties of the Hamiltonian 
under investigation. This determines the parameters Pi 
and Cij. An interaction length z is chosen according to 
the desired emulation time t. 

To translate the theoretical design into a real exper- 
imental system, the geometry of the circuit is analysed 
in a numerical electromagnetic design suite, in our case 
RSoft [35]. This program uses finite-difference algo- 
rithms to find solutions to Maxwell's equations in di- 
electrics and can be used to optimise a set of physical 
parameters (core diameter and refractive index contrast) 
which will approximate the desired evolution. This pro- 
vides the link between the Hamiltonian evolution of the 
quantum system under study, and the physical properties 
of the experimental system. 

We now describe the two-photon evolution in the opti- 
cal system. Two temporally-indistinguishable input pho- 
tons 1^) = aja^ |0), for i ^ j, give the joint detection 
probability ^ of finding the state b\bj |0) in output 
modes k and / [16]: 

ri,, = {blbihkh = ^^\Ui,kU,,i + U,M\ (4) 

1 + Ok,l 

The matrix = {F^^} then describes the two-photon 
probability distributions in all combinations of output 
modes. In contrast, two distinguishable photons will 
evolve independently and obey the statistics of Bernoulli 
trials. The corresponding output probability distribution 
F^ I takes the form, 

= TZJ- d^^.fcf^^-.'l' + \Wj,k?) ■ (5) 

In the case of photons and electro-magnetic fields, F^ ^ 
represents an intensity correlation matrix F^ ^ = {Ikh)- 
The components F^ ^ and F^ ^ will differ by the factor 
2Yle[{Ui^kUj^iY{Ui^iUj^k)\i which encompasses the quan- 
tum nature of indistinguishable walkers. 

To quantify the overlap fidelity between 
two probability distributions, we use S = 




Mueller matrix 

The Mueller-matrix array M describes the transfer 
properties of the circuit in terms of the Stokes param- 
eters <S, which describe the polarisation state of an elec- 
tromagnetic field [29]. The component j is the 4x4 
Mueller matrix describing the coupling from input waveg- 
uide j to output waveguide i. From the corresponding 
output Stokes parameter Sij for the input state 1^)^, 
the real- valued parameter \Uij\^ was calculated as the 
output transmission component in the \H). subspace. 

Table 1 shows the calculated Mueller-matrix array M 
of the quantum walk circuit. The matrix not only al- 
lows us to determine the evolution matrix U in the \H) 
subspace, but also quantifies the polarisation-dependent 
coupling and birefringence observed in the circuit, as seen 
in figure [4] 
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